For orbit inside Earth's 
magnetosphere, start here: 



CREME96 ROAD 

For orbits at or beyond Auxiliary Programs you may 
geosynchronous, start here: wish to run: 
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Input Particle Rigidity in GV, 
arrival directions (i.e. zenith 
angles 6 min &9 m&x and 
azimuth angles (p min & <p max ). 
number of arrival directions per 
orbital step (M^,), and 
geomagnetic activity level (Kp 
and Dst). 



L 



Choose number of steps per 
each orbit (N sp0 ) of the Earth 
(e.g. N spo = 200 steps), set 
orbit step counter N = 0 and 
number transmitted N tr = 0, 
and specify required accuracy 
of numerical integration 
(typically 10 s ). 
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Input orbital parameters: 
inclination, apogee, perigee, 
initial time (to), orbital duration 
(T), initial displacement of 
perigee from ascending mode, 
initial longitude of ascending 
node, and initial displacement 
from the ascending mode. 
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Calculate orbital period (P) 
from input orbital parameters 
using standard orbit 
generation routines (e.g. 
Adams et al., 1986). 
Calculate: 



^^N^M^T/P 



and St = P/N. 
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Find spacecraft's orbital location 
for time t= to +N<5t using 
numerical solution of standard 
equations for orbit generator, e.g. 
those given on pp. 43 - 48 of 
(Sterne, 1960), as implemented 
in the CREME orbit generator 
(Adams etal., 1986). Set particle 
trajectory time { = t and .time 
step 8 tf=0. Set M = 0. 



Differential equation for integrating is 
the Lorenz equation in magnetic field, 
using "backward" solution (Q replaced 
by -Q, r 7 replaced by -r 7 ) of: 

F = m ydv/dt 7 = QvxB 

B = B, GRF (r t {) + B Tsyg89 (Kp,r, {) 
+ 5B cxlcndedTsyg g 9 (Kp,Dst, r/) 

NOTE: In the "backward" solution, 
the actual particle's final position is 
the initial position along the numerical 
integration of the "backward" particle 
trajectory, and the "backward" 
particle's initial velocity is in the 
opposite direction as the actual 
particle's final velocity. 

See text for description of B-field 
calculations. 
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Set r = particle's final position and v 
= particle's final velocity vector. 
These are determined from the 
spacecraft's latitude, longitude, 
and altitude, combined with the 
particle's rigidity, final 0, and 
final <p. 



Generate particle's final 6 and 
<j£> for each orbital location, using 
random number techniques, i.e. 
X! and x 2 are two distinct outputs 
from a random number generator. 

COS0=COS0 max + 

x, (COS0 min -COS0 max ) 

<P=9 min+X 2 ((j9 max - (pmin) 
M = M+ 1 

NOTE: many modern computers 
have built in random number 
generator functions. For 
computers which do not, one can 
use the subroutines for random 
number generation from (Press et 
al., 1992). 



Evaluate B-field and d v/dr 7 at location 
r and time \! along the particle's 
trajectory, i.e. t 7 is the time along the 
particle's trajectory. 
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Specifically, perform Bulirsch-Stoer 
numerical integration, using formulae & 
algorithms as presented on pp. 718-725 
of (Press et al., 1992). The outputs are 
the particle's radial position r, velocity 
v, and time t 7 =t+<5 1 7 after the numerical 
integration step; and an estimate of <5t / 
for the next particle-trajectory 
integration time-step. 

NOTE: 8 \! is re-calculated within each 
execution of the adaptive Bulirsch-Stoer 
integration. 



Evaluate if particle trajectory 
encountered solid Earth 
(forbidden), entered atmosphere 
(forbidden), crossed 
magnetospheric boundary 
(transmitted), or still within 
magnetosphere. 
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If transmitted, N ff = N„ + 1 
If forbidden, Nn, = !*«,+ 1 



If within magnetosphere. 

If transmitted or forbidden, 
N = N+ 1 
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If M < M.^ 
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IfN<N mM 
If N > N mK 
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Calculate Geomagnetic 
Transmission (GT) for this 
rigidity (R) and output the result: 

GT(R) = N ff /N m „ 

NOTE: N^N^-N,, 



^^ r: y47W,f t/u^&ML /jZ?j 
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CREME96 Main Me 
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See brief instructions at the bottom of this page. See How to Run CREME96 for more details. 
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1 . To CREATE a User Request Ffo : 

o Click on the Edit button by the corresponding routine, leaving the User Request File blank. (If necessary, use 

the blank space in the pulldown menu to clear a name already typed there.) 
o Click on the "GO" bar. 
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To EDIT an* existing User Request File : 
o Click on the Edit button by the corresponding routine, 
o Select a User Request File from the appropriate pull-down menu, 
o Click on the "GO" bar. 



To RUN one or more routines: 
o If an Edit button has been clicked "on", use the RESET bar to clear it. 
o Select User Request File(s) from the appropriate pull-down menu(s). 
o Click on the RUN button(s) in that line(s) 
o Click on the "GO" bar. 



4. To access CREME96 Utilities: 

o Click on the RUN button on that line 
o Click on the "GO" bar. 

Note: You can create or edit only^HifeJJserRequest File at a time. But you can submit several files for running 
(sequentially, in the order shown in thetaLneTby^hecking more than one "RUN" button and then clicking on the "GO" 
bar. 

Questions? Comments? Send mail to tylka@crs2.nrl.navv.mil . 
You may terminate your CREME96 session at any time. 
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